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What is claimed is: 

1 . An apparatus for computing multiple integral of a multidimensional integrand 
function A to be integrated with using a vector map f with unbounded support which 
converts an m (m>l) -dimensional vector having real number components into an m- 

5 dimensional vector having real number components, by which a multidimensional density 
function p for the limiting distribution resulting from repeatedly applying the map f to an 
m-dimensional vector u is analytically solvable, said apparatus comprising: 

a first storage unit which stores an m-dimensional vector x = (x l5 x 2 , . . . , xj; 
a second storage unit which stores a scalar value w; 
1 0 a first computing unit which computes a vector x 1 = f(x) = (x^, x' 2 , . . . , x'J resulting 
from applying said vector map f to said vector x being stored in said first storage unit; 

a second computing unit which computes a scalar value w r = A(x)/p(x) based on 
said vector x being stored in said first storage unit and said scalar value w being stored in 
said second storage unit; 
1 5 an update unit which updates the value stored in said first storage unit by storing 
said vector x' computed by said first computing unit on said first storage unit, and updates 
the value in said second storage unit by adding said scalar value w' computed by said 
second computing unit to a value to be stored in said second storage unit; and 

an output unit which computes a scalar value s = w/(c+l) based on said scalar value 
20 w being stored in said second storage unit when the number of update times by said 
update unit becomes c (c>l), and outputs said scalar value s as a result of the multiple 
integral. 

2. The apparatus according to claim 1, wherein said scalar value stored in said 
second storage unit first is a result from dividing a value resulting from applying said 

25 function A to said m-dimensional vector stored in said first storage unit first by a value 
resulting from applying said density function p to said m-dimensional vector stored in 
said first storage unit first. 
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3 . The apparatus according to claim 1, wherein said scalar value stored in said 
second storage unit first is 0, and 

said output unit computes a scalar value s* = w/c, and outputs said scalar value s' as 
the result of the multiple integral instead of said scalar value s. 
5 4. The apparatus according to claim 1 further comprising: 

a convergence rate obtainer which obtains convergence rate of scalar values 
sequentially output by said output unit while varying said number of update times c for 
each of plural vector maps g l9 g 2 , . . . , g k (k>2) which are prepared as said vector map f; 
a vector map selector which refers to the convergence rates obtained by said 
1 0 convergence rate obtainer, and selects a vector map g h (l<h<k) which shows fastest 
convergence rate; and 

an output controller which controls said output unit to output said scalar value with 
using said vector map g h as said vector map f and the number of update times c' (c'>c) 
instead of said number of update times c. 
15 5. The apparatus according to claim 1, wherein a multidimensional density 
function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u)))), . . . 
resulting from applying said vector map f to a predetermined m-dimensional vector u = 
(u l9 u 2 , . . . , uj for equal to or greater than 0 times, satisfies the following property of: 

20 p(u) = n=o m Pi<A); 

p,(iO ~ M~ a+a} for u, -» -co; 
~ a, |uj- (1+a) for u, — +co; 
(a>0, l<i<m, c_>0, c +1 >0) 
6. The apparatus according to claim 5, wherein said vector map f is defined as 
25 f(u) = (fiOiJ, f 2 (u s ), . . . , 4(10) 

by a function f;(t) = g^d, t)/d 1 (d >0) which is defined in l<i<m, and said map g, is 
defined by any one of the following maps cpj (l<j<8) and a natural number (n>2), as 
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follows: 

g 1 (( Pj (0)) = (p j (n i 0); 
^(0) = -sgn(tan0) / |tan0| 1/a ; 
(p 2 (9) = -sgn(tan0) x |tan0| 1/a ; 
5 q> 3 (8) = -sgn(cos0) / |tan0| 1/a ; 
q> 4 (0) = -sgn(cos0) x |tan0| 1/a ; 
q> 5 (0) = sgn(cos0) / [tan0| 1/a ; 
q> 6 (8) = sgn(cos0) x |tan0| 1/a ; 
(f) 7 (0) = sgn(sin0) / |tan0| 1/a ; 
10 cp 8 (9) = sgn(sin0) x |tan0| 1/a ; 
sgn(t) = 1 for t > 0; 
sgn(t) = 0 fort = 0; 
sgn(t) = -l fort<0 

7. The apparatus according to claim 5 further comprising: 
1 5 a convergence rate obtainer which defines said map f for each of plural positive 
numbers q u . . . , % (k>2) prepared as an invariable a, and obtains convergence rates of 
the scalar values sequentially output by said output unit while varying said number of 
update times c; 

a positive number selector which refers the convergence rates obtained by said 
20 convergence rate obtainer, and selects an integer q h (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said map f with using said positive number as 
said invariable a, and controls said output unit to output said scalar values with using the 
number of update times c' (c'>c) instead of said number of update times c. 
25 8. The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map g, with using plural ones of said 
maps <p-, and obtains convergence rates of the scalar values sequentially output by said 
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output unit while varying said number of update times c; 

a map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects one of said maps cpj which shows the fastest convergence rate; 
and 

5 an output controller which defines said map & with using said map q>j selected by 
said map selector, and controls said output unit to output the scalar values with using the 
number of update times c' (c'>c) instead of said number of update times c. 

9. - The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map g, relating to each of plural 
1 0 natural numbers p 1; p 2 , . . . , p k (k>2) as said natural numbers n 13 and obtains convergence 
rates of the scalar values sequentially output by said output unit while varying said 
number of update times c; 

a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number p h (l<h<k) which shows the 
1 5 fastest convergence rate; and 

an output controller which defines said natural number map & with using said 
natural number p h as said natural number % and controls said output unit to output the 
scalar values with using the number of update times c' (c'>c) instead of said number of 
update times c. 

20 10. The apparatus according to claim 1, wherein said output unit computes said 
scalar value s each time said update unit carries out update, compares said latest scalar 
value s with the former scalar value which is computed at former update, and outputs said 
latest scalar value s if a result of the comparison satisfies a predetermined condition for 
terminating the computation. 

25 11. A method for computing multiple integral of a multidimensional integrand 
function A to be integrated with using a vector map f with unbounded support which 
converts an m (m>l) -dimensional vector having real number components into an m- 
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dimensional vector having real number components, by which a multidimensional density 
function p for the limiting distribution resulting from repeatedly applying the map f to an 
m-dimensional vector u is analytically solvable, 

a first storage unit which stores an m-dimensional vector x = (x l5 x 2 , . . . , xj, and 
5 a second storage unit which stores a scalar value w, 

said method comprising the steps of: 

computing a vector x' = f(x) = (x' l5 x' 2 , . . . , x'^ resulting from applying said vector 
map f to said vector x being stored in said first storage unit; 

computing a scalar value w' = A(x)/p(x) based on said vector x being stored in said 
1 0 first storage unit and said scalar value w being stored in said second storage unit; 

updating the value stored in said first storage unit by storing said vector x' computed 
by said first computing unit on said first storage unit, and updating the value in said 
second storage unit by adding said scalar value w' computed by said second computing 
unit to a value to be stored on said second storage unit; and 
1 5 computing a scalar value s = w/(c+l) based on said scalar value w being stored in 
said second storage unit when the number of update times by said update unit becomes c 
(c>l), and outputting said scalar value s as a result of the multiple integral. 

12. The method according to claim 11, wherein said scalar value stored in said 
second storage unit first is a result from dividing a value resulting from applying said 

20 function A to said m-dimensional vector stored in said first storage unit first by a value 
resulting from applying said density function p to said m-dimensional vector stored in 
said first storage unit first. 

13. The method according to claim 11, wherein said scalar value stored in said 
second storage unit first is 0, and 

25 said outputting step computes a scalar value s' = w/c, and outputs said scalar value s' 
as the result of the multiple integral instead of said scalar value s. 

14. The method according to claim 1 1 further comprising the steps of: 
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obtaining convergence rate of scalar values sequentially output by said outputting 
step while varying said number of update times c for each of plural vector maps &, 
■ ■ ■ > §k (k>2) which are prepared as said vector map f; and 
referring to the convergence rates obtained by said convergence rate obtainer, and 
5 selecting a vector map g h (l<h<k) which shows fastest convergence rate, and 

said outputting step outputs said scalar value with using said vector map g h as said 
vector map f and the number of update times c 1 (c'>c) instead of said number of update 
times c. 

15. The method according to claim 1 1 , wherein a multidimensional density 
1 0 function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u))), . . . 
resulting from applying said vector map f to a predetermined m-dimensional vector u = 
(u„ u 2 , . . . , uj for equal t0 01 gr eater man 0 times > satisfies the following property of: 
p(u) = FU m PXuD; 
15 fcOO ~ cjuj -0 "* 0 for U;-^ -co; 

pXu,) ~ c_, |uj- (1+a) for u, — +go; 
(a>0, l<i<m, c_>0, c + >0) 

16. The method according to claim 15, wherein said vector map f is defined as 

f(u) = ftCuO, f 2 (u s ), . . . , fJyO) 
20 by a function f^t) = g^d, t)/d t (d>0) which is defined in l<i<m, and said map g, is 
defined by any one of the following maps y } (l<j<8) and a natural number n, (n>2), as 
follows: 

g 1 (cp j (e)) = cp j (n i e); 
9^6) = -sgn(tan0) / |tan9| 1/a ; 
25 cp 2 (G) = -sgn(tan9) x |tan0j I/a ; 

q> 3 (6) = -sgn(cos0) / ltan6| 1/a ; 
(p 4 (0) = -sgn(cose) x |tan0| 1/a ; 
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(p s (0) = sgn(cos0)/|tan6| 1/a ; 
cp 6 (0) = sgn(cos0) * |tanG| 1/a ; 
(p 7 (6) = sgn(sin0) / |tan6| 1/a ; 
q> 8 (0) = sgn(sin0) * |tan0| 1/a ; 
5 sgn(t) = l fort>0; 
sgn(t) = 0 fort = 0; 
sgn(t) = -1 for t < 0 

17. The method according to claim 15 further comprising the step of: 

defining said map f for each of plural positive numbers q u q 2 , . . . , % (k>2) prepared 
1 0 as an invariable a, and obtaining convergence rates of the scalar values sequentially 
output by said output unit while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a positive number 
(l<h<k) which shows the fastest convergence rate; and 

defining said map f with using said positive number q^ as said invariable a, and 
1 5 controlling output of said scalar values with using the number of update times c' (c'>c) 
instead of said number of update times c. 

18. The method according to claim 16 further comprising the steps of: 
defining said map g, with using plural ones of said maps q> p and obtaining 

convergence rates of the scalar values sequentially output by said output step while 
20 varying said number of update times c; 

referring to the obtained convergence rates, and selecting one of said maps q>j which 
shows the fastest convergence rate; and 

defining said map g, with using said selected map cpj selected, and controlling output 
of the scalar values with using the number of update times c' (c'>c) instead of said number 
25 of update times c. 

19. The method according to claim 16 further comprising the steps of: 
defining said map g t relating to each of plural natural numbers p„ p 2 , ■ - ■ , Pk ( k ^ 2 ) 
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as said natural numbers % and obtaining convergence rates of the scalar values 
sequentially output by said outputting step while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a natural number p h 
(l<h<k) which shows the fastest convergence rate; and 
5 defining said natural number map g, with using said natural number p h as said 
natural number n 15 and controlling output of the scalar values with using the number of 
update times c' (c'>c) instead of said number of update times c. 

20. The method according to claim 1 1 , wherein said outputting step computes 
said scalar value s each time said updating step carries out update, compares said latest 

1 0 scalar value s with the former scalar value which is computed at former update, and 
outputs said latest scalar value s if a result of the comparison satisfies a predetermined 
condition for terminating the computation. 

21. A computer readable recording medium storing a program for computing 
multiple integral of a multidimensional integrand function A to be integrated with using a 

1 5 vector map f with unbounded support which converts an m (m>l) -dimensional vector 
having real number components into an m-dimensional vector having real number 
components, by which a multidimensional density function p for the limiting distribution 
resulting from repeatedly applying the map f to an m-dimensional vector u is analytically 
solvable, said program causes a computer to function as: 

20 a first storage unit which stores an m-dimensional vector x = (x l3 x 2 , . . . , xj; 
a second storage unit which stores a scalar value w; 

a first computing unit which computes a vector x' = f(x) = (x\, x' 2 , . . . , x' J resulting 
from applying said vector map f to said vector x being stored in said first storage unit; 

a second computing unit which computes a scalar value w' = A(x)/p(x) based on 
25 said vector x being stored in said first storage unit and said scalar value w being stored in 
said second storage unit; 

an update unit which updates the value stored in said first storage unit by storing 
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said vector x' computed by said first computing unit on said first storage unit, and updates 
the value in said second storage unit by adding said scalar value w' computed by said 
second computing unit to a value to be stored on said second storage unit; and 

an output unit which computes a scalar value s = w/(c+l) based on said scalar value 
5 w being stored in said second storage unit when the number of update times by said 
update unit becomes c (c>l), and outputs said scalar value s as a result of the multiple 
integral. 

22. The recording medium according to claim 2 1 , wherein said scalar value 
stored in said second storage unit first is a result from dividing a value resulting from 

1 0 applying said function A to said m-dimensional vector stored in said first storage unit first 
by a value resulting from applying said density function p to said m-dimensional vector 
stored in said first storage unit first. 

23 . The recording medium according to claim 21 , wherein said scalar value 
stored in said second storage unit first is 0, and 

1 5 said output unit computes a scalar value s' = w/c, and outputs said scalar value s' as 
the result of the multiple integral instead of said scalar value s. 

24. The recording medium according to claim 2 1 , wherein said program further 
causes said computer to function as: 

a convergence rate obtainer which obtains convergence rate of scalar values 
20 sequentially output by said output unit while varying said number of update times c for 
each of plural vector maps g l7 (k>2) which are prepared as said vector map f; 

a vector map selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a vector map g h (l<h<k) which shows fastest 
convergence rate; and 

25 an output controller which controls said output unit to output said scalar value with 
using said vector map g h as said vector map f and the number of update times c' (c'>c) 
instead of said number of update times c. 
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25. The recording medium according to claim 21, wherein a multidimensional 
density function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u))), . . . 
resulting from applying said vector map f to a predetermined m-dimensional vector u = 
5 (u l3 u 2 , . . . , uj for equal to or greater than 0 times, satisfies the following property of: 

p(u) = n 1 =o m P i (u i ); 

p^U;) ~ c_, |uj~ cl+a) for u, — * -co; 
p^uj ~ c_, M -0 "* 0 for Ui —»• +oo; 
(a>0, l<i<m, c_>0, c + >0) 
10 26. The recording medium according to claim 25, wherein said vector map f is 
defined as 

f(u) = (f 1 (u 1 ),f 2 (u s ),...,f m W) 
by a function f,(t) = g,(d, t)/d, (d>0) which is defined in l<i<m, and said map g t is 
defined by any one of the following maps cpj (l<j<8) and a natural number n, (n>2), as 
15 follows: 



g 1 (9 J (e)) = ( Pj (n 1 6); 


<Pi(e) = 


-sgn(tanG) / |tanG| 1/a ; 


<fe(e) = 


-sgn(tanG) x |tan9| 1/a ; 


(P 3 (0) = 


-sgn(cosG) / |tan0| 1/a ; 


<?M = 


-sgn(cosG) x |tan0| 1/a ; 


q> 3 (e) = 


sgn(cosG) / |tanG| 1/a ; 


q> 6 (e) = 


sgn(cosG) x |tanG| 1/a ; 


(p 7 (9) = 


sgn(sinG) / |tanG| 1/a ; 


q> 8 (e) = 


sgn(sinG) x |tanG| 1/a ; 


sgn(t) = 


= 1 for t > 0; 


sgn(t) = 


= 0 for t = 0; 


sgn(t) = 


= -1 for t < 0 
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27. The recording medium according to claim 25, wherein said program further 
causes said computer to function as: 

a convergence rate obtainer which defines said map f for each of plural positive 
numbers q l5 q^ . . . , % (k>2) prepared as an invariable a, and obtains convergence rates of 
5 the scalar values sequentially output by said output unit while varying said number of 
update times c; 

an integer selector which refers the convergence rates obtained by said convergence 
rate obtainer, and selects a positive number (l<h<k) which shows the fastest 
convergence rate; and 

1 0 an output controller which defines said map f with using said positive number q^ as 
said invariable a, and controls said output unit to output said scalar values with using the 
number of update times c' (c'>c) instead of said number of update times c. 

28. The recording medium according to claim 26, wherein said program further 
causes said computer to function as: 

1 5 a convergence rate obtainer which defines said map g, with using plural ones of said 
maps (p j; and obtains convergence rates of the scalar values sequentially output by said 
output unit while varying said number of update times c; 

a map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects one of said maps q>j which shows the fastest convergence rate; 

20 and 

an output controller which defines said map g, with using said map ^ selected by 
said map selector, and controls said output unit to output the scalar values with using the 
number of update times c' (c'>c) instead of said number of update times c. 

29. The recording medium according to claim 26, wherein said program further 
25 causes said computer to function as: 

a convergence rate obtainer which defines said map g, relating to each of plural 
natural numbers p l3 p 2 , . . . , p k (k>2) as said natural numbers n i; and obtains convergence 
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rates of the scalar values sequentially output by said output unit while varying said 
number of update times c; 

a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number p h (l<h<k) which shows the 
5 fastest convergence rate; and 

an output controller which defines said natural number map g t with using said 
natural number p h as said natural number n^ and controls said output unit to output the 
scalar values with using the number of update times c 1 (c'>c) instead of said number of 
update times c. 

1 0 30. The recording medium according to claim 21, wherein said output unit 

computes said scalar value s each time said update unit carries out update, compares said 
latest scalar value s with the former scalar value which is computed at former update, and 
outputs said latest scalar value s if a result of the comparison satisfies a predetermined 
condition for terminating the computation. 

15 31. The recording medium according to claim 2 1 , wherein said recording medium 
is a compact disc, a floppy disk, a hard disk, a magneto-optical disk, a digital versatile 
disc, a magnetic tape, or a semiconductor memory. 



